Fermi surfaces and quantum oscillations in underdoped high- Tc superconductors 

YBaaCugOe.s and YBa2Cu408 
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We study underdoped high-Tc superconductors YBa2Cu306.5 and YBa2Cu408 using first- 
principles pseudopotential methods with additional Coulomb interactions at the Cu atoms, and 
obtain Fermi-surface pocket areas in close agreement with measured Shubnikov-de Haas and de 
Haas-van Alphen oscillations. With antiferromagnetic order in Cu02 planes, stable in the calcu- 
lations, small hole pockets are formed near the so-called Fermi-arc positions in the Brillouin zone 
which reproduce the low- frequency oscillations. A large electron pocket, necessary for the negative 
Hall coefficient, is also formed in YBa2Cu306.5, giving rise to the high-frequency oscillations as 
well. Effective masses and specific heats are also calculated and compared with measurements. Our 
results highlight the important role of magnetic order in the electronic structure of underdoped 
high-Tc superconductors. 

PACS numbers: 71.18.-(-y, 74.25.Jb, 74.72.-h, 74.25. Ha 



The normal-state electronic structures of the under- 
doped high-Tc superconductors have been studied for 
more than twenty years, but the Fermi-surface (FS) 
topology is still only partially understood [l|, 0, 0, 0, H, 
i, 0, i i, E [il IH El Q E H El- An impor- 
tant observation is the disconnected FSs l|, Q , namely 
Fermi arcs, observed in angle-resolved photoemission 
spectroscopy (ARPES), which initiated intense investi- 
gations about whether the FSs are really disconnected 
arcs or closed pockets of which one side is hardly vis- 
ible in ARPES. Recently, in contrast to having Fermi 
arcs, de Haas-van Alphen (dHvA) oscillations in the 
magnetization and Shubnikov-de Haas (SdH) oscillations 
in the resistance^, 0, S S S Hi observed in ortho-H 
YBa2Cu306.5 and YBa2Cu408 suggest well-defined close 
pockets in the FS of underdoped high-Tc cuprates. The 
measured oscillations for YBa2Cu306.5 are a dominant 
one at 500±20 T with a satelhte at 1650±40 T [H*!, and 
more recently a dominant oscillation at 540±15 T with 
satellites at 450±15 T, 630±40 T, and 1130±20 T H. For 
YBa2Cu408, oscillation at 660±15 T is observed [T, 

The measured dHvA and SdH oscillations provide ex- 
treme cross-sectional areas of closed FS pockets [isj, but 
they alone are not enough to identify the shapes and 
locations of the pockets. Thus, a quantitative theoreti- 
cal calculation of the FS geometry can be useful to de- 
termine the FS topology. First-principles calculations 
based on the density functional theory (DFT) approach 
have been performed for YBa2Cu306.5 and YBa2Cu408 
[1; [H U I J but the calculated FSs could not explain the 
oscillation measurements. In contrast to the meansure- 
ments, reported DFT calculations predict only FS pock- 
ets much larger than 500 T, and do not obtain the 
electron-type carriers implied by the observed negative 
Hall coefficients 121. 



According to model calculations [H 0, [3 
antiferromagnetic (AFM) order or a d-density wave may 
result in small pockets in regions of the FS. Although us- 
ing DFT calculations, one may examine static magnetic 
order using spin-density functional theorvfioj: as yet, no 
magnetic order has been considered in the reported DFT 
calculations of YBa2Cu306.5 and YBa2Cu408. 

In this Letter, we present, for the first time, 
first-principles spin-density functional calculations of 
YBa2Cu306.5 and YBa2Cu408 with a Coulomb repul- 
sion U at Cu sites which yield FSs consistent with the 
dHvA and SdH measurements. It is shown that, with 
physically reasonable U values, the AFM order in the 
Cu02 planes reconstructs the FS and produces pock- 
ets with sizes consistent with the measured frequencies. 
Moreover, the calculated FS of YBa2Cu306.5 has a large 
electron pocket which explains the observed negative Hall 
coefficients. In addition, cyclotron effective masses and 
specific heats are calculated and compared with exper- 
iments. Our results support the possible importance of 
magnetic order in the electronic structures of underdoped 
high-Tc cuprates. 

Our present work is based on ab-initio pseudopoten- 
tial density-functional calculations with pseudo-atomic 
orbitals to expand the electronic wavefunctions [20| . 
Coulomb interaction at Cu d orbitals, parameterized by 
U and J 21, [23|, is added to the local (spin) density 



approximation [L(S)DA-|-[/]. With experimental atomic 
structures [2K2JI, we minimize the total energy with re- 
spect to the magnetic moments of Cu atoms in the Cu02 
planes and CuO chains to consider the possibility of AFM 
order. Our results are that YBa2Cu306 is an AFM in- 
sulator and YBa2Cu307 is a non-magnetic metal. 

Using the LDA+J7 method with no magnetic order, we 
obtain the electronic structure for YBa2Cu306.5 (Figs, la 




FIG. 1: (Color online) Electronic structures in ortho-II 
YBa2Cu306.5 with Coulomb interaction at Cu d orbitals {U 
— 6.0 eV and J = 1.0 eV). (a) The band structure and (b) 
the FS obtained by the LDA+C/ method with no magnetic 
order, (c) The band structure and (d) the FS obtained by 
the LSDA+(7 method with a checker-board-type AFM order 
in the Cu02 planes. In (b) and (d), the FSs are drawn in the 
Brillouin zone of a real-space unit cell (0.383 x 0.387 nm^) 
containing a Cu atom on each Cu02 plane. Blue (red) lines 
are hole pockets and open orbits (electron pockets), (e) FS 
pockets in (d). The a, a' , and /3 pocket areas are 485, 621, 
and 1450 T, respectively, (f) Comparison with experimental 
FS pocket areas (630 T Q and 1650 T ;|l) in dashed lines. A 
FS area of 1 nm"'^ corresponds to a frequency of 105 T [l^ . 

and lb), which are in good agreement with previous 
calculations!^ [13]. The FS has only large hole pockets 
and open orbits (Fig. lb); however, this is not in agree- 
ment with the observed quantum oscillations. 

When magnetic order is considered in the LSDA-|-[/ 
calculations for YBa2Cu306.5, a checker-board-type 
AFM order is stabilized in the Cu02 planes, and this 
drastically changes the electronic structures (Figs. Ic and 
Id). With U ^ 6.0 eV and J = 1.0 eV, each Cu atom 
in the Cu02 plane has 0.48 Bohr magneton (/is), while 
the CuO chains are still non-magnetic. The FS (Fig. Id) 
now consists of small hole pockets (a and a' indicating 
the two largest ones) and a large electron pocket {(3). The 
calculated pocket areas, which are not very sensitive to 
U and J around the used values, are in good agreement 
with the experimental observations (Fig. If). This shows 
that the AFM order may be a way to quantitatively ex- 
plain the measured dHvA and SdH frequencies. 

Figure Id shows that the hole pockets {a and a') are 
located at (±^, ±^), close to the positions of the Fermi 



FIG. 2: (Color online) FS pocket sizes versus fcz and cyclotron 
effective masses in YBa2Cu306.5. (a) FS pocket areas perpen- 
dicular to the fez-axis when the Fermi level is at the charge- 
neutrality level, (b) Hole and (c) electron pocket areas and 
their cyclotron effective masses as functions of the Fermi-level 
shift (AEf) from the charge-neutrality level. The indexes a, 
a' , and f3 are the same as in Fig. 2d. A FS area of 1 nm"'^ 
corresponds to a frequency of 105 T [l^ . 

arcs m the ARPES dataP, 01 This supports the idea 
that slow AFM fluctuation'25j may form hole pockets 
near (±^, ±-|), with their shapes possibly modified to 
form the arcs because of short-range fluctuation (isl. Il3| . 

Figure Id also shows that the electron pocket (/?) is 
much more anisotropic than the hole pockets (a and a'). 
This arises because the electron pocket is derived from 
the CuO chain and the Cu02 plane, while the hole pock- 
ets come from the Cu02 plane only. To have an isotropic 
resistivity as observed in experiments at temperature be- 
low 80 K [2^ , we find that the electron mean free path 
along the chain direction should be about one quarter of 
that perpendicular to the chain because of directional 
difference in the group velocity. 

For more detailed comparison. Fig. 2a shows the kz 
dependence of FS pocket areas. The three largest ex- 
treme areas of hole pockets are 485 and 621 T at kz—Q 
and 708 T at kz—n/c (Fig. 2a), which overestimate by 
about 15 % the experimental low frequencies, 450, 540, 
and 630 T Q . The largest extreme area of the electron 
pocket (/3) is 1450 T at kz=0 (Fig. 2a), and it underesti- 
mates by 12 % the experimental high frequency, 1650 T 
Q. If the Fermi level is shifted to higher energy, hole 
pockets would shrink and the electron pocket would ex- 
pand (Figs. 2b and 2c). With a Fermi-level shift (AEp) 
of 4 meV above the charge-neutrality level, the extreme 
pocket sizes become 441, 564, and 652 T for the holes and 
1519 T for the electron, respectively, resulting in closer 
agreement with experimental results. 

We calculate the cyclotron effective masses (Figs. 2b 
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FIG. 3: (Color online) Hall coefficients for YBa2Cu306.5. (a) 
Hall coefficient versus magnetic field (B) at temperature T 
= 1.5 K, with and without the SdH oscillation (red and blue 
lines, respectively), (b) Hall coefficient versus temperature at 
B = 60 T without considering the SdH oscillation. 



and 2c) and the electronic contribution to the normal- 
state specific heat from the LSDA+t/ electronic struc- 
ture of YBa2Cu306.5. The obtained cyclotron effective 
masses are 1.78 times the free electron mass (me) for the 
a' pocket and 1.88 me for the (3 pocket. These values 
are smaller than measured values, 1.78 ~ 1.9 rUe for the 
low frequency and 3.8 me for the high frequency [1, 0, 0], 
but they are consistent with experiments in the sense 
that the effective mass of the low-frequency oscillation 
(from the a' pocket in our result) is smaller than that of 
the high-frequency oscillation (from the (3 pocket in our 
result). The calculated Sommerfeld coefficient for the 
normal-state specific heat is 9.28 mJ mol~^ which 
slightly underestimates the experimental value of 10 niJ 



mol""'^ [27[. The differences between our values and 
the measured ones may originate from many-body effects. 

The presence of the electron pocket in our FS (Fig. Id) 
definitely opens a chance of a negative Hall coefficient, 
but it alone is not sufficient since the total numbers of 
holes and electrons in our calculation are equal to each 
other to represent a charge-neutral stoichiometric sam- 
ple. Since the Hall coefficient is inversely proportional 
to the net charge of the carriers, a slight imbalance of 
the two types of carriers would result in a relatively large 
Hall coefficient. With AE'f = 4 meV, as discussed above 
for the oscillation frequencies, we can obtain a negative 
Hall coefficient of -25 mm^C"! at 70 T (Fig. 3ay^close 
to the experimental value of about -30 mm 

Figure 3 shows the Hall coefficients obtained by semi- 
classical transport theory within the relaxation-time ap- 
proximation, assuming temperature-dependent but field- 
independent mean free paths (A) and A£'i? — 4 meV. At 
50 K, we assume isotropic A's for holes, which are 40 nm 
for the a and a' pockets and 10 nm for the other smaller 
pockets, and anisotropic A for electrons (/3), which are 
20 and 80 nm for motion along and perpendicular to the 
CuO chain, respectively. At 1.5 K they are increased 
to ten times the values at 50 K. These values of A's are 
adjusted to show a theoretical reproduction of the exper- 
imental data although they are quite a bit larger than 



those estimated from oscillation amplitudes. With the 
assumed A's, the calculated Hall coefficient is negative at 
high magnetic field and low temperature (Fig. 3a), be- 
coming positive at high temperature (Fig. 3b), as in the 
experiment (12i] . 

The SdH oscillations are displayed in Fig. 3a by modi- 
fying the conductivity tensor <Tij to include effects of the 



Landau levels 28| . The above mentioned mean free paths 
(A) are used for Oij itself; however, a reduction of A's by 
a factor of 0.05 is assumed for the modification factor of 
Oij for quantum oscillations, yielding SdH oscillation am- 
plitudes close to experiments[l2]. This may suggest that 
the charge carriers in the material lose their quantum co- 
herence much faster than their classical linear momenta, 
but it is beyond the scope of our present work to justify 
the assumed A's. In our results, the low-frequency oscil- 
lations (from the hole pockets) are much stronger than 
the high-frequency oscillation (from the electron pocket) 
since the average mean free time is larger for the holes 
than for the electrons even with the assumed A's because 
of the difference in their group velocities. Thus, the dom- 
inant oscillation in the Hall coefficient (Fig. 3a) originates 
from the hole pockets although the Hall coefficient itself 
is negative at high field due to the electron pocket. The 
calculated SdH oscillations grow with the magnetic field 
(Fig. 3a), as observed experimentally |,12i| . 

Compared with model calculations considering mag- 
netic ffuctuations 13j, [ij] , our results show that the pres- 
ence of CuO chains in YBa2Cu306.5 is important for ex- 
plaining the high-frequency quantum oscillation and the 
negative Hall coefficient. Since our results are based on a 
static long-range magnetic order stable in the LSDA-|-?7 
method, fluctuations in real materials may modify the 
FS. As discussed above, one possibility is the evolution 
of the small-size hole pockets (a and a' in Fig. Id) to 
arcs, as proposed by the model calculations. 

While the d-density-wave theory predicts hole pock- 
ets larger than electron pockets, our result predicts an 
electron pocket larger than hole pockets, and assigns FS 
pockets to the observed frequencies oppositely. Thus, in 
our work, the observed major frequency originates from 
small-size hole pockets while the negative Hall coefficient 
is due to a large-size electron pocket. 

For YBa2Cu408, as in the case of YBa2Cu306.5, we 
obtain very different FSs with and without AFM order 
in the Cu02 planes (Fig. 4). With the LDA+U method 
with U = 3.1 eV and J = 0.8 eV for all Cu d orbitals, 
we obtain a FS consisting of large hole pockets and open 
orbits (Fig. 4b) which is consistent with previous first- 
principles calculations (2^. When the magnetic order is 
considered by the LSDA-|-?7 method with the same U 
and J, AFM order is stable in the Cu02 planes with 
0.22 But each Cu atom, and the FS consists of small 
hole pockets (a and a' indicating the two largest ones), 
open orbits, and small electron pockets (/3), as shown 
in Fig. 4d. The calculated FS pocket areas are 721 T 
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FIG. 4: (Color online) Electronic structures in YBa2Cu408 
with Coulomb interaction at Cu d orbitals {U — 3.1 eV and J 
= 0.8 eV). (a) The LDA+!7 band structure and (b) the FS in 
the non-magnetic case, (c) The LSDA+C/ band structure and 
(d) the FS with a checker-board-type AFM order in the Cu02 
planes. In (b) and (d), the FSs are drawn in the momentum 
space corresponding to a real-space unit cell (0.384 x 0.387 
nm'^) containing a Cu atom on each Cu02 plane. Blue (red) 
lines are hole pockets and open orbits (electron pockets), (e) 
FS pockets in (d). The a, a' , and /3 pocket areas are 721, 729, 
and 685 T, respectively, (f) Comparison with the measured 
FS pocket area (660 T 0,0) i^i the dashed line. 



(a), 729 T {a'), and 685 T (/?), which are close to the 
measured value 660 T 0, [§] i overestimating it by about 
10 % or less. Contrary to YBa2Cu306.5, the calculated 
FS pocket sizes in YBa2Cu408 are sensitive to U and J 
around the used values. Calculated cyclotron effective 
masses, 0.45 TOg for holes and 0.52 me for electrons, are 
much smaller than the measured values of 2.7 3.0 me 0, 

, but calculated Sommerfeld coefficient for the normal- 
state specific heat, 6.97 mJ mol~^ K~^, is rather close to 
the experimental value of 9 mJ mol^^K"^ (soj . 

In summary, we have studied the electronic struc- 
tures of YBa2Cu306.5 and YBa2Cu408 by the LSDA-I-C/ 
method, and the results yield FS topologies fully consis- 
tent with quantum oscillation measurements. It is shown 
that the magnetic order in the Cu02 planes may ex- 
plain quantitatively the dHvA and SdH oscillation fre- 
quencies, the negative Hall coefficients, and the specific 
heat. These results point to the importance of magnetic 
order for understanding the electronic structures of the 
underdoped high- Tc cuprates. 
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